Robustness of the noise-induced phase synchronization 
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Abstract 

We show that a wide class of uncoupled limit cycle oscillators can be in-phase synchronized 
\^ \ by common weak additive noise. An expression of the Lyapunov exponent is analytically derived 
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chronization can be achieved in a broad class of oscillators with little constraint on their intrinsic 
property. On the other hand, the leaky integrate-and-fire neuron oscillators do not belong to this 
^ ' class, generating intermittent phase slips according to a power low distribution of their intervals. 
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Populations of nonlinear oscillators can be found in a variety of phenomena, including 
laser array [iL semiconductors |2], chemical reactions [3|, society of living organisms j^, 
and neurons S] . In many of these systems, the phases of oscillations can precisely coincide 
owing to mutual interactions among oscillators. Alternatively, a strong periodic input may 
synchronize independent oscillators through the entrainment to the common input. In either 
cases, external or internal noise sources may disturb the phase synchronization, and therefore 
have long been considered to exert a negative influence on the precise temporal relationship 
between oscillators. 

This view, however, has been challenged recently jo, l2|. Pikovsky studied, in his pio- 
neering work, a population of circle maps stimulated by impulse inputs at discrete random 
times, and found that the common noise can induce stable phase synchronization. Since 
the noise-driven synchronization does not depend on the intrinsic frequency of oscillators, 
it differs from the entrainment to an external periodic input. Evidence is accumulating for 
the common noise-induced synchronization in several biological and physical systems. For 
instance, an ensemble of independent neuronal oscillators may be synchronized by a fluc- 
tuating input applied commonjy to all of them. This is suggested by experimental studies 
of neural information coding 8|, in which the reproducibility of spike firing was tested for 
a repeated application of the same fluctuating or constant input current. Interestingly, the 
reproducibility of the output spike trains was much higher for the fluctuating input than 
for the constant one ja Isl lll| , indicating a high temporal precision of the spike responses 
to noisy input. In ecological systems, common environmental fluctuations such as climate 
changes may synchronize different populations separated by a large geographical distance 



ll[. In fluid dynamics, a common turbulent flow may generate a synchronized motion of 



floating particles 
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All these findings indicate an active role of noise in synchronization of non-interacting 
dynamical elements. It remains, however, unclear whether the noise induced phase synchro- 
nization is specific to a limited class of oscillators, or can be generalized to a broad class 
of oscillators. In this study of a general class of limit-cycle oscillators, we show that com- 
mon additive noise, even if it is weak, can induce phase synchronization regardless of their 
intrinsic properties and the initial conditions. Using the phase reduction method which is 
applicable to an arbitrary oscillator [3|, we analytically calculate the Lyapunov exponent of 
the synchronizing state and prove that the exponent is non-positive as long as the phase- 



dependent sensitivity is different iable up to the second order. In addition, we investigate 
the scahng laws that appear in the dynamics of the relative phase when the perfect phase 
synchronization is deteriorated by a discontinuous phase-dependent sensitivity or oscillator- 
specific noise sources. 

Population of A^ identical nonlinear oscillators driven by common additive noise are de- 
scribed as 

Xi = F{X,) + m, (1) 

where i = 1,- ■ ■ ,N and $,{t) is a vector of Gaussian white noise. The elements of the vector 
are normalized as < ^i(t) >= and < ^i(t)^m(s) >= 2DimS(t — s), where D = {Dim) is a 
variance matrix of the noise components. Because all the oscillators are identical and do not 
interact with one another, we can study the phase synchronization of the entire population 
in a reduced system of two oscillators. Regarding the common noise as a weak perturbation 
to the deterministic oscillators, the phase reduction method gives the following dynamical 
equations of the phases: 

4)i = uo + Z{<Pi) ■ ^, (2) 

where u is an intrinsic frequency of the unperturbed oscillators. Z is the phase-dependent 
sensitivity defined as Z{(j)) = gradx4>\x=Xoi^)^ where Xo(0) is the unperturbed limit cycle 
solution determined by F{X). We assume that Z is differentiable at least to the second 
order, although Z can be discontinuous for such oscillators that have discontinuous periodic 
solutions (e.g., integrate-and-fire neurons). As we will see later, the discontinuity of Z can 
significantly affect the noise-driven synchronization. To ensure the validity of the phase 
reduction, the weak noise must satisfy the condition uj ^ \Dim\- 

Equation (J2I) implies that the synchronizing solution described as 0i(t) = 02 (^) is absorb- 
ing, i.e., once two oscillators synchronize, they always remain synchronizing. Since the area 
of the phase space is limited (0 < 0i,02 < Svr), the phase variables starting from arbitrary 
initial phases can reach a neighborhood of the synchronizing solution with a finite probability 
in a finite time. To prove that the synchronizing solution is stable against perturbations, we 
analytically calculate the Lyapunov exponent A of the solution. We note that the stochastic 
equation Q should be interpreted as a Stratonovich differential equation, since the phase 
reduction method assumes the conventional variable translations in differential equations. 
To remove the correlation between and ^, we translate Eq. (j2)) into an equivalent Ito 



differential equation [13[: 

(i>^ = u; + Z'{<f),)DZ{<j)i) + Z{<f),) ■ ^, (3) 

where dash denotes differentiation with respect to 0. Suppose that the two phases have an 
infinitesimally small difference tp = (f)2 ~ (j^i- Then, linearization of Eq. Q with respect to 
ip gives 

^ = {{z\cp)Dzm' + z\(P) ■ m. (4) 

where obeys Eq. Q. By introducing a new variable y = log{ip), Eq. (JD) is further rewritten 
as 

y = {Z'DZy - {Z'DZ') + Z' ■$,. (5) 

Since the Lyapunov exponent is defined as YwiT^ooiyiT) — y{0))/T, the long time average 
of right hand side of Eq. (j3)) coincides with A. Replacing the long time average with the 
ensemble average with respect to ^, we can represent A as 

A = < {Z'DZy - {Z'DZ') + Z' ■i>^ 
= < {Z'DZy - {Z'DZ') >g 
= / P,tm{Z'DZy-{Z'DZ')]dcp. (6) 

Here, the second line follows from < Z'{(f)) ■$, >^= 0, which holds in Ito stochastic processes. 
Pst is a steady distribution function of (p described as 

Psti<P) = ^ / ^ ^ ^^P (^(^) - yi<t>))d^^ (7) 

with an effective potential V"(0) = — J ^ ziXoziX ^-"^ ^^^ ^ normalization constant C. 
Fortunately under the assumption of weak noise, u ^ \Dim\, Pgt is reduced to a constant 
function Pst{,(p) = l/(27r). By substituting this into the last line of Eq. ©, and noting that 
the first term vanishes due to the periodicity of Z, we finally obtain the following formula: 

A = -— / Z'DZ'dcp < 0, (8) 

271" Jo 

where the equality holds if Z' is a constant function. Since the variance matrix is positive 

definite, A is non-positive. This implies that the phase synchronization induced by common 

noise is stable in an arbitrary oscillator system regardless of the detailed oscillatory dynamics, 

as long as Z is differentiable. 
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FIG. 1: The common noise- induced phase synchronization, (a) The phase difference between 
two oscillators shows an exponential decay ffuctuating around the theoretical behavior (solid line) 
derived from Eq.Q. Different curves correspond to different realizations of the random driving 
force ^(t). Here, Z((f)) = sin0, uo = \ and D = 0.1. (b) The Lyapunov exponent of the above 
synchronizing solution is shown as a function of the common-noise intensity (circles) . The solid line 
shows an analytical result, (c) The Lyapunov exponent is calculated for the synchronizing solution 
of the Stuart-Landau oscillator. Numerical results are shown for both original (filled circles) and 
reduced (empty circles) oscillator systems. The solid line represents a theoretical result. The 
parameter values are set as cq = 2, C2 = 1, Dn = D22 = D and D12 = D21 = 0. 

To confirm the above analytical results, we numerically solve Eq. (J2I) to obtain the Lya- 
punov exponent for a specific choice of Z, i.e., Z{(f)) = sm{(f)). The phase difference between 
the two oscillators driven by common additive noise shows an exponential decay, and the 
decay constant well agrees with the analytical result (Fig. ^). Consistent with equation 
(jHl), the magnitude of the negative Lyapunov exponent increases in proportion to the noise 
intensity D (Fig. lb). In order to confirm the validity of the phase reduction method, we 
employ the Stuart-Landau oscillators and compare the Lyapunov exponent derived from 
Eq. (jHJ with that calculated numerically in the original oscillator system (Fig. Ic). For the 
Stuart-Landau oscillator described as A = (1 + ico)A — (1 + 2C2)|v4py4, the phase sensitiv- 
ity can be explicitly given as Z{(f)) = (— C2cosc<j0 — sina;0, — C2sina;0 + cosuj(f))/uj, where 
a; = Co - C2 J2|. 
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FIG. 2: The intermittent phase shps induced by uncommon additive noise sources to the oscillators 
defined by w = 1 and Z((j)) = sini;^. The noise intensities D = 0.1 and d = 0.001. (a) The 
distribution of the phase difference over a sufficiently long time, (b) The time evolution of the phase 
difference exhibits intermittent phase slips, (c) The distribution of the inter-phase-slip intervals 
shows a power-low decay with an exponent of —3/2 (fitted by a solid line) and an exponential 
cutoff resulted from the additive noise term in Eq. (fTUI) . 

In practical situations, the individual oscillators may be influenced by additional, 
oscillator-specific noise sources. They may appear from the fluctuations intrinsic in the 
oscillators, or a lack of perfect coincidences in the common driving noise. To discuss the 
influences of additional noise, Eq. (P) is modified to 



X, = F{Xi)+^it) + ri,it), 



(9) 



where the uncommon noise sources r]i are normalized as < Vi,iit) >= ^^^ < 
Vi,i{'t)Vj,m{s) >= 2dimSij6{t — s). Linearization of Eq. (jH)) gives the stochastic equation 
of the phase difference as 



^ 



{Z'{D + d)Zy + Z'-{^ + 



^ + Z ■ (772 - r7i) 



(10) 



Since both multiplicative and additive factors fluctuate, Eq. ()1()|1 is regarded as a multiplica- 
tive stochastic process with additive noise, which has been studied in variety of fields 14| . 
The steady distribution function of Eq. (fTIH) exhibits a power-low decay in a middle range of 
ip (Fig. 2a), implying that for most of the time the phases stay near phase synchronization 
even in the presence of uncommon noise sources. Simulations of Eq. (jU)) reveal that the 
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FIG. 3: The unstable phase synchronization of integrate-and-fire models stimulated by a common 
noise source. The phase-dependent sensitivity of this oscillator is discontinuous, and the noise- 
driven synchronization is not ensured. 1 = 2 and no uncommon noise, d = 0. (a) The Lyapunov 
exponent A is plotted as a function of the common noise strength D. (b) The phase difference shows 
significant jitters, diffusing due to intermittent phase slips. Here and below, D = 0.1. (c) The time 
evolution of the two phasess displays temporarily synchronizing {t < 170) and desynchronizing 
{t > 170) states, (d) The distribution of the inter-phase-slip intervals shows an exponential decay 
with an exponent of —3/2 (fitted by a solid line). 



system is trapped in the phase synchronizing state for certain intervals between intermittent 
phase slips (Fig. 2b). It is known that the intermittent bursts are characteristic to the 
stochastic processes driven simultaneously by multiplicative and additive noise sources |l5[, 
and that the intervals between neighboring bursts obey a power low distribution with an 
exponent of —3/2. In Fig. 2c, the inter-slip intervals of the present phase dynamics obeys 
a power low distribution of the same exponent. 

So far, we have assumed that the phase-dependent sensitivity Z is a. continuous function 
of the phase. However, the phase reduction method does not ensure the continuity of Z, and 
some oscillators do not have this property. For example, an integrate-and-fire neuron oscilla- 



tor, which is described hj i) = I—v with a renewal condition v{t) = 1 ^ hm^^o+ v(t+T) = 0, 
is frequently used for modeling neuronal activity, but it has the following discontinuous Z: 

Z(0) = yexp(^), O<0<27r, (11) 

where u = 27r/(log/ — log (/ — 1)). As shown in Fig. 3a, numerical integrations of Eq. Q for 
the Z given in Eq. (jllll show positive Lyapunov exponents, with the magnitudes increased 
with the intensity of the common noise. The phase difference does not decay exponentially, 
but fluctuates around the synchronizing state between the intermittent phase slips (Fig. 3b). 
Figure 3c displays the synchronized time evolution of the phase variables that is terminated 
by an abrupt phase slip at t ~ 170: After that, the two phases are desynchronized until 
they recover the phase synchronization (not shown). As in the previous case, the inter-slip 
intervals obey a —3/2 power- low distribution (Fig. 3d). This intermittency is essentially the 
same as the on-off intermittency of chaotic oscillators just before the onset of synchroniza- 
tion, thus associated with positive Lyapunov exponents |la|. Note that the discontinuity 
of the phase sensitivity and the resultant positive Lyapunov exponents are inherent in the 
leaky integrate-and-fire model. For example, the Hodgkin-Huxley model has a continuous 
and different iable phase sensitivity, thus yielding a negative Lyapunov exponent and a stable 
phase synchronization in response to common noise (results not shown). 

We briefly argue the relationships between the present study and two previous studies. 
In the stochastic resonance, the ability of an excitable system in detecting a weak signal 
can be optimized by noise of suitable intensity [1/3] • The present study also argues the role 
of noise in improving the response reliability. However, here the improvement is achieved 
by the precise temporal coincidences between oscillators, whereas the stochastic resonance 
only enhances the response probability without caring the exact timing of events. Thus, the 
stability or Lyapunov exponent is not a central issue in the stochastic resonance, and the 
two studies deal with qualitatively different phenomena. Second, some chaotic osciUators 
exhibited phase synchronization when they were driven by common additive noise 
However, it remained unknown whether this type of synchronization may appear in a broad 
class of, either chaotic or non-chaotic, oscillator systems. In this paper, we have proven that 
such a synchronization can be induced in a wide class of limit-cycle oscillators. A unified 
treatment of limit cycle oscillators and chaotic oscillators is awaited, as they may share 
many characteristic properties of the common noise-induced synchronization. 
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In conclusion, independent limit cycle oscillators can be synchronized by weak, common 
additive noise regardless of the detailed oscillatory dynamics and the initial phase distri- 
butions. The stability of this synchronizing solution only requires the presence of a second 
derivative of the phase-dependent sensitivity, so the solution can exist in a broad class of 
oscillators. The leaky integrate-and-fire oscillators do not belong to this class of oscillators, 
and show no perfect synchronization. The analytical evaluation of the Lyapunov exponent 
remains open for further studies for oscillators possessing non-differentiable phase sensitivity. 
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